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We present the diagrammatic technique for calculating the free energy of the Hermitian one- 
matrix model to all orders of 1 /N expansion in the case where the limiting eigenvalue distribution 
spans arbitrary (but fixed) number of disjoint intervals (curves). 



1 Introduction 



^1- 

■ Matrix models and their so-called multisupport (multicut) solutions became again important recently 
due to studies in M = l 1 SUSY gauge theories by Cachazo, Intrilligator and Vafa [1], [2] followed by 
the proposal of Dijkgraaf and Vafa [3] to calculate the low energy superpotentials using the partition 

"^T* 1 function of multicut solutions. These solutions, being known already for a long time [4, 5], got a new 

Oh. insight in [6]. 

The leading order of the 1/iV-expansion in the matrix size of the matrix model is described by 
. the semiclassical tau-function of the so-called universal Whitham hierarchy [7] (see the details about 
^ ' one-matrix and two- matrix cases in [8] and [9]). 

One may also consider more general solutions to matrix models, identifying them with generic 
solutions to the loop (Schwinger-Dyson, or Virasoro) equations [10], to be the Ward identities satisfied 
by matrix integrals [11]. 

We consider solutions to the loop equations admitting multi-matrix integral representations [6, 12, 
13]. These solutions are associated with families of Riemann surfaces and form a sort of a basis in the 
space of all solutions to the loop equations [13] (like the finite-gap solutions form a similar basis in the 
space of all solutions to an integrable hierarchy). They can be distinguished by their "isomonodromic" 
properties — switching on higher matrix model couplings and 1 /TV-corrections does not change the 
family of Riemann surfaces, but just reparameterizes the moduli as functions of these couplings. On 
the side of integrable hierarchies, these solutions must satisfy some equations of multicomponent 
hierarchies [14], not just the Toda chain hierarchy equations, which are satisfied only by a one-cut 
solution. 

Much progress has been made on the way to constructing explicit solutions of matrix models in 
the large- N expansion since [15]. It was shown (primarily for the one-cut case [16]) that the variables 
convenient for describing solutions in the genus expansion are branching points of the corresponding 
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Riemann surface and the so-called moments of the potential. These variables, being expressed highly 
nonlinearly in times and occupation numbers, allow reducing the problem to solving some algebraic 
relations. Using this technique, the two- matrix model was solved in the large- N limit in the one- 
cut case [17], this technique was extended to the multicut case [9], [18] , and the subleading term of 
1/iV-expansion was found in [19]. 

In [20] , the consistent diagrammatic technique for finding the loop means (multiresolvents) in the 
1/iV-expansion of the Hermitian one-matrix model (1MM) was proposed. The main result of the 
present paper is the generalization of this technique to describe the free energy contributions of all 
orders in 1/N for multicut solutions of the 1MM. 

In Sec. 2, we describe the general properties of multi-cut solutions of 1MM and the solution of the 
loop equation in the large- iV limit. 

In Sec. 3, we first present necessary notions of the Riemann geometry describing the Bergmann 
bidifferentials and the loop insertion operator (which acts from the space of s-differentials to the space 
of (s + l)-differentials) and then describe the diagrammatic technique of [20] in a more convenient 
terms than in the original presentation. 

In Sec. 4, we introduce the operator H, which reduces the degree of forms by one and is in a 
sense inverse to the loop insertion operator. This eventually makes possible to present the genus k 
contribution to the free energy through the action of H on the one-point resolvent (for all k except 
k = 1; for the latter case, the answer has been already obtained, see [21], [19]). 



2 Matrix models in 1/iV-expansion 

2.1 Loop equation and resolvents 

Let us consider the formal integral of the 1MM [22] 

/ DXe-^ V{X) =e-^ (1) 
Jn 



NxN 



where V(X) = J2n>i tnX n and fi = to/N is a formal expansion parameter. The integration in (1) goes 
over NxN Hermitian matrices, and for generic potentials, the integration may go over curves in the 
complex plane of each of N proper variables, which are eigenvalues of X. The topological expansion 
of the Feynman diagrams series is then equivalent to the expansion in even powers of U for 

oo 

F = F{h,to,t u t2,...) = Y J tf 9 ~ 2 F g - (2) 

9=0 

Customarily, to = is the scaled number of eigenvalues. It is convenient, but not compulsory, 
to assume the potential V(p) to be a polynomial of the fixed degree m + 1. In general, to have an 
algebraic curve we need V'(p) to be a rational function [23] 

The averages corresponding to partition function (1) are defined as usual: 

(F(X)) = 1 J n ^DXF(X) exp (-±trV(*)) , (3) 
and it is convenient to use their formal generating functionals: the one-point resolvent 

as well as the s-point resolvents (s > 2) 

w <»- ■■■?.}=* 2 - t <tr y + v • tr i t : > r- = (tr -b - *r ) < 5 > 

fcl,...,fc s =l Pi ---Ps \ Pi A p s VL/ conn 
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where the subscript "conn" pertains to the connected part. 

These resolvents are obtained from the free energy T through the action 

d d dT 



w(pi,...,p s ) = -h 



dv( Ps ) dv( Ps ^) dv( Pl ) " 

d d d 

W( P1 ), (6) 



of the loop insertion operator 2 



av(p,) av{p,. 1 ) dvfa 



dv( P ) - at 
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Therefore, if one knows exactly the one-point resolvent for arbitrary potential, all multi-point resolvents 
can be calculated by induction. In the above normalization, the genus expansion has the form 

oo 

W( Pl ,...,p s ) = Y,fr 2g W g (pi,...,p s ), s>l, (8) 

9=0 

which is analogous to genus expansion (2). 

The first in the chain of the loop equations of the 1MM is [10] 

duJ V '^W{uj) = W{pf + h 2 W(p,p). (9) 



Cj, 2m p — uj 

Here and hereafter, Ct> is a contour encircling clockwise all singular points of W(oj), but not the 
point uj = p; this contour integration acts as the projection operator extracting the coefficient of the 
term Using Eq. (6), one can express the second term in the r.h.s. of loop equation (9) through 
W(p), and Eq. (9) becomes an equation on one-point resolvent (4). 

Substituting genus expansion (8) in Eq. (9), we find that W g (p) for g > 1 satisfy the equation 

9~ 1 pi 

(k-2W (p)) W g (p) = J2 W g ,{p)W g „ g ,{p) + ^-jW fl _i(p), (10) 
where K is the linear integral operator 

In Eq. (10), W g (p) is expressed through only the W gi (p) for which g, L < g. This fact permits developing 
the iterative procedure. 

The solution W±(p) to the loop equation in the multicut case was first found by Akemann [24]. 3 He 
also managed to integrate it in order to obtain the free energy T\ in the two-cut case. The genus-one 
partition function in the generic multi-cut case was proposed in [25, 26], where it was observed that 
the Akemann formula coincides with the correlator of twist fields (that produce cuts on the complex 
plane and give rise to the hyperelliptic Riemann surface as a cover of complex plane) computed by 
Al.Zamolodchikov [27], with some corrections due to including star operators (introduced by G.Moore 
in [28]) into consideration. In [21, 19], the genus-one correction was derived by solving the loop 
equation, which generalize the Akemann result for partition function to arbitrary number of cuts. 

2 Although this operator contains all partial derivatives w.r.t. the variables tk's, below we introduce additional variables 
Si and, therefore, we use the partial derivative notation here. 

3 The universal critical behavior of the corresponding correlation functions was discussed in [5]. 
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2.2 Solution in genus zero 



The form of loop equation (9) is based exclusively on the reparameterization invariance of the matrix 
integral, which retains independently on the details of eigenvalue density distribution. In the 1MM 
case at N — ► oo, the eigenvalues fill in some segments in complex plane, dependently on the shape of 
potential V(X). For polynomial potentials, the number of segments is finite and the contour Cv of 
integration in (10) encircles a finite number n of disjoint intervals 

n 

v = U^ai-i'A^i], Mi < M2 < • • • < M2„. (12) 
i=i 

Recall that as 



W g {p)\v^oo = t j5 gfi + O{l/ P % (13) 



all W g (p) must be total derivatives, 



W »< p ) = WM?«> 9>0. (14) 



Inserting genus expansion (8) into loop equation (9), we obtain 

d\ V'(X) 



Cj, 2iri p — A 



Wo(\) = (W (p)Y (15) 



for genus zero and (10) for higher genera. Given Wq(j>), one can then determine W g (jp) for g > 1 
iteratively genus by genus provided the operator (K — 2Wo(p)) (see (11)) can be inverted uniquely. 

One can solve Eq. (15) for the planar one-point resolvent Wq(p) as follows. Deforming the contour 
in Eq. (15) to infinity, we obtain 



(W (p)) 2 = V'(p)W (p) + I t^^Wa). (16) 

JCoo 2tti p — a 

The last term in the r.h.s. is 



a degree-m — 1 polynomial to be fixed later, and the solution to (16) is then 

Wo(p) = \V(jp) - yV'( P y + 4P m ^(p) = lv» - yip), (18) 

where the minus sign is chosen in order to fulfill asymptotic condition (13) and the function y(p) is 
defined as follows. For the polynomial potential of power m + 1, the resolvent Wq(p) is a function on 
complex plane with m cuts, or on a hyperelliptic curve y 2 = V'{p) 2 + AP m ^i(p) of genus g = m — 1. 
For generic potential V(X) with m —>■ oo, this curve may have an infinite genus, but we can still 
consider solutions with a finite, fixed number n of cuts filled by eigenvalues. For this, we separate the 
smooth part of the curve introducing 

2n 



yip) = M(p)y(p), and y 2 (p) = J[ a=1 {p ~ Ma) (19) 



with all branching points jj, a distinct. The variable y defines therefore the new, reduced Riemann 
surface, which plays a fundamental role in our construction. In what follows, we still assume M(p) to 
be a polynomial of degree m — n, keeping in mind that n is always finite and fixed, while m > n can 
be chosen arbitrarily large. 
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From now on, we distinguish between images of the infinity at two sheets — physical and unphysical — 
of hyperelliptic Riemann surface (19) respectively denoting them 00+ and oo_. By convention, we set 
y\ P ^oo + ~ P n , and M(p) is then 4 

1 V'(w) 

M(p) = --re Soo+ dw- \^r- (20) 

2 (w — p)y(w) 

Inserting this solution in Eq. (18) and deforming the contour back, we obtain the planar one-point 
resolvent with an n-cut structure, 

Let us now discuss the parameter counting. Rewriting (16) as quadratic equation for Wo(p) we 
had to imply that polynomial (17) depends on some unknown parameters to be fixed later. Indeed, 
the dependence on coefficients of (17) can be retranslated into the dependence on filling fractions 

*-h&« x >-L£i M ^ (22) 

where Ai, i = l,...,n — 1 is the basis of A-cycles on the hyperelliptic Riemann surface (19) (we 
may conveniently choose them to be the first n — 1 cuts). Adding the (normalized) total number of 
eigenvalues 

h = j c y( A )^~ = res °°+ y( x ) dX ( 23 ) 

to the set of Si, we obtain n parameters, which exactly matches the number of coefficients of (17) in 
the nondegenerate case where m = n. We assume the occupation number for the last, nth cut to be 
to — Y%=i $i = S n . 5 If n < m, we consider variables (22) as conditions to which we add, following 
(13), the asymptotic conditions 

1 / dX X k V'(X) , , 

which provide n + 1 equations for 2n constants u Q . 

In the planar limit of matrix models, filling fractions (22) can be considered as independent pa- 
rameters of the theory, since the jumps between different cuts are suppressed as non-perturbative 
corrections in h. (Obviously, no parameters Si arise in the one-cut case.) In particular, this imposes 
restrictions 

ww) s ' = ' • ' - wtf°-°- (25) 

In accordance with the original matrix model formulation (see e.g. [4], [24]) one must look for the 
genuine minimum of the matrix model effective action with respect to all variables: Si together with 
the times tj. This implies vanishing of partial derivatives of J-q with respect to Si, and one can find 
(see [4], [3], and [29]) that these derivatives are differences of chemical potential on disjoint cuts, equal 
to the integrals over dual S-cycles on (19): 



dSi 



I y{X)dX = ^i (26) 

JBi 



Note that the geometric definition of is not modular invariant, i.e., it depends on the choice of A- 
and .B-cycle basis on (19). Under the change of homology basis, transforms in accordance with the 
duality transformations [30]. The higher-genus corrections become also scheme dependent: choosing 
Si or Hi as independent variables, we obtain different expressions for the genus-one free energy. 

4 By a standard convention, res^dx/x — —1, and the direction of the integration contour at the infinity point therefore 
coincides with the direction of contour for integrals over Ct> and over the set of ^4-cycles, see below. 

J It is sometimes convenient to consider S n instead of to as a canonical variable. However, in all instants we use S n , 
we specially indicate it for not confusing S n with the "genuine" filling fraction variables Si, i = 1, . . . ,n — 1. 



5 



3 Calculating resolvents. Diagram technique 



In this section, we rederive the diagrammatic technique of [20] in a more concise form of differentials 
on Riemann surface (19). We choose Si to be independent variables. Our main goal is to invert loop 
equation (10) to obtain the expression for Wk(p) for any k > 1. 

3.1 Bergmann bidifferential and two-point resolvent 

Let us consider the Bergmann kernel (canonically normalized bidifferential in Fay's terminology) which 
is the bi-differential on a Riemann surface T, g that is symmetrical in its arguments F,Q £ S 9 and has 
the only singularity at the coinciding arguments where it has the behavior (see [31], [32]) 

B(P, Q) = ( (r(p) _ 1 r(Q))2 + \S B (P) + c(l)) dr(P)dr(Q), (27) 

in some local coordinate t(P); Sb(P) is the Bergmann projective connection and we fix the normal- 
ization (the possibility to add symmetrical bilinear forms composed from holomorphic 1-differentials 
in variables Q and P) claiming vanishing all the integrals over ^-cycles of B(P, Q): 

<f B(P,Q) = 0, for i = l,...,g. (28) 

We then have the following standard Rauch variational formulas relating B(P,Q) with other objects 
on a (general, not necessarily hyperelliptic) Riemann surface: 

° B(P,Q) = ]-B(P,\ Ma ])B(\fjL a ],Q), (29) 



<9/i 
and 

/ B(P,Q) = 2-iridwi(P), (30) 

JBi 

where fi a is any simple branching point of the complex structure, then, by definition, in the vicinity 
of n a , 

B(P,Q)\ Q ^ a =B{PM) ( ^L_ + 0(^7^)dg) , (31) 
and dWi(P) are canonically normalized holomorphic differentials: 

/ dwi(P) = Sij. (32) 
JAj 

We also introduce the 1-form dEQ t Q (P), which is the primitive of B(P,Q): 

[Q dr(P) 
dE Q>Qo (P)= J B(P,Z), dE Q>Qo (P)\ P ^ Q = { _ r j (Q) + finite . (33) 

Then, obviously, 

/ dE Q>Qo (P)=0. (34) 
JAi 

The form dEQ t Q Q (P) is single- valued w.r.t. P on the Riemann surface and multiple-valued w.r.t. the 
variable Q: from (30), 

dE ( P )Q+§ B ,Q^ P ) = 2iridwi(P)+dE QtQo (P), 
dE{P) Q+§ l Qo {P) = dE Q , Qo {P). 
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The reference point Qo does not play any role in the construction; we keep it only for the consistency 
of notation. 

These quantities are related to the Prime form E(P, Q): 

B(P,Q) = d P d Q logE(P,Q), dE QiQo (P) = ^ log ^®\ , (35) 

E{P, Qo) 

where the Prime form is defined in the standard way. Let us consider the Jacobian J, which is a 
^-dimensional torus related to the curve £ 9 . Recall that the Abel map S 9 ^J: Q — > x = j Jq q duji j , 
where Qo is a reference point, set into the correspondence to each point Q of the complex curve 
the vector in the Jacobian, and we also introduce the theta function 6[ a ](x) of an odd characteristic 
[a] that becomes zero at x = 0. We introduce the normalizing functions (1/2-differentials) h a {x) 
determined for the subset of x G J that are image points of Q G S 9 

= E^(o)<mq). 

i=l 1 

The explicit expression for the Prime form E(x, y) that has a single zero on the Riemann surface E g 
then reads 

K{x)h a (y) 

We can now express the 2-point resolvent Wo(p, q) in terms of B(P, Q). Let us denote p and p the 
points on the respective physical and unphysical sheets. Then, 

8V " W =- B (P, 5 )- B (P,5)- ** 



since it has double poles with unit quadratic residues at p = q and p = q. The 2-point resolvent (18) 
is nonsingular at coinciding points; therefore, 

^ = -l( B <P>9)- B <P>*))> W 

and 

W (p,q) = -B(p,q). (38) 

3.2 The iterative procedure 

We can determine higher genus contributions iteratively by inverting genus expanded loop equation 
(10). All multi-point resolvents of the same genus can be obtained from W g (p) merely by applying 
the loop insertion operator • 

Claiming all the higher free energy terms F g to depend only on /j, a and a finite number of the 

moments M& \ which are derivatives of (k — l)th order of the polynomial M(p) at branching points, 
allows no freedom of adding the terms depending only on to and Si to F g . 

We first prove that the operator K — 2Wo{p) acting on 

c/ E ^ M ^) m - <39) 

where f{q) is the combination of resolvents in the r.h.s. of (10) whose only singularities are poles of 
finite orders at \i a and the point p lies outside Ct>, just gives f(p). 

For this, we first note that as V'(£) has no square root singularities at fi a , instead of evaluating 
the integral over Cx>, we can just evaluate residues on Riemann surface (19) w.r.t. the variables £ and 
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q. (This effectively corresponds to doubling the integral over Cp and it gives additional factor 1/2.) 
We therefore denote the contour that encircles all branching points of y(£) and assume that the 
contour C(«) of integration w.r.t. the variable q lies inside . 

We now implement the Riemann bilinear identities. The contributions of A- and -B-cycle integra- 
tions vanish due to (25) and (30). Following (18), we substitute 2(W (0 + y(0) for V'(£) in (11) 
and calculate the part with Wo(£) taking the contour of integration to the infinity; the integrand with 
dEq,q (0 1S regular, and the only contribution comes from poles at £ = p and £ = p; this contribution 
exactly cancels the term —2Wq(p). 

The remaining integral is 

_r d^d^vi^m , , f(g) 

7 P >C(0> C ( 9 ) 2ni2nip-tm q > qo ^> 2y(q)- 

Since y{0/y(0 = is regular at £ = fi a , pushing the integration contour through we 

obtain zero upon the integration over £, and the only nonzero contribution comes from the simple pole 
of dE qqo {£) ~ l/(£ — q) at £ = q. Terms with y{q) in the numerator and denominator then cancel each 
other and the remaining integral w.r.t. q can be taken pushing the contour again to infinity evaluating 
the residues at q = p and q = p, which gives \{f(p) + f{p)) = f{p) for a function analytic in p outside 
the branching points. 

This provides a basis for the diagrammatic representation for resolvents in 1MM [20]. Let us 
represent the form dE qM) (p) as the vector directed from p to q, the three-point vertex as the dot 
in which we assume the integration over q, • = § ^ 2y{q) ' an< ^ ^ ne Bergmann 2-form B(p,q) as a 
nonarrowed edge connecting points p and q. The graphic representation for a solution of (10) then 
looks as follows. Representing the multiresolvent W g >(pi, . . . ,p s ) as the block with s external legs and 
with the index </, we obtain [20] 




(40) 

which provides the diagrammatic representation for Wk{pi, . . . ,p s ). The multiresolvent Wk(pi, . . . ,p s ) 
can be presented as a finite sum of all possible connected graphs with k loops and s external legs and 
with only three- valent internal vertices (the total number of edges is then 2s + 3k — 3, and we assume 
s > 1 for k > 1 and s > 3 for k = 0). We also set arrows on some (exactly 2k + s — 2) of edges for 
the arrowed edges to constitute the rooted tree subgraph with all arrows directed from the root. That 
means that we choose one of the external legs, say, p\ (the choice is arbitrary due to the symmetry 
of Wk(pi, ■ ■ ■ ,p s )), to be the root vertex the tree starts with; for each three- valent vertex there must 
exist exactly one incoming edge of the tree subgraph, all external edges (except the root edge) are 
nonarrowed, and all internal nonarrowed edges either start and terminate at the same vertex (we then 
associate with such an edge the Bergmann differential B(p,p)) or connect two different vertices such 
that there exists the directed path composed by arrowed edges connecting these two vertices. At each 
internal vertex we have the integration § C ( q ) ^ > while the arrangement of the integration contours 
at different vertices is prescribed by the arrowed subtree: the closer is a vertex to the root, the more 
outer is the integration contour. 

We now demonstrate the consistency of this diagram technique by calculating the action of loop 
insertion operator (7) on its elements. 

We first calculate the action of d/dV(r) on B(P, Q). Using (29), we represent this action through 
the action of partial derivatives in \x a subsequently calculating the latter from relation (37). Let 

y{x)dx\ x ^ tla = y{\n a ])y/x - ^ a dx + 0(^x - n a ) 3 dx. 
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Then, since 



we have 



and, therefore 



dy(p)dp 



dV(r) 



1 dp dfi c 

-oJ/(M- 



d^g _ 2B([n a ],q) 
dV(r) ~ y(\p a ]), 



(41) 



9 R/pnl _f P(P> a ])P([^Q)P(K],r) 

#v(r) ^ ~~ ~1 v(M • (42J 

For expressing this in terms of the differentials, we integrate one of the Bergmann bidifferentials in 
order to obtain the 1-differential d£ , Q i Q (^). 6 Note that the local variable in the vicinity of \i a is 
= y/x — /j, a , this gives the additional factor 1/2, and we eventually obtain 

9 B(PO)-^T res B(P,C(x))dE^ )tQo (Q)B(C(x),r) 
From this relation it obviously follows that 



(n , ^ dE P , Q Mx))dE axl Q (.Q)B(Z(x),r) 
dE P>Qo (Q) = ^ res^ o„(^M^c M 



dV(r) ^ "« 2y(xi(x))di(x) 

and the last quantity to evaluate is 

d 1 B(p,r) 



(44) 



dV(r)2y(p) 2y\p) 



(45) 



Note that the point P in (44) is outside the integration contour. (This is irrelevant in formula (43) 
as the Bergmann bidifferential has zero residue at the coinciding arguments, in contrast to dEpQ (Q). 
Multiplying by l/(2y(P)) both sides of (44), using (45), and pushing the integration contour through 
the point P, we observe that the contribution of the simple pole of dEp t Q (£(x)) at = P cancels 
exactly the variation of l/(2y(P)). We therefore attain the prescribed contour ordering and can 
graphically present the action of d/dV(r) as 



Wir) Q — *P = Q^+P, dW) Q P = <2— ^ — P = Q 



(46) 

In the second case, it is our choice on which of edges to set the arrow. Recall, however, that the points 
P and Q were already ordered as prescribed by the diagram technique. That is, if "P — > Q" , we must 
choose the first variant and if "Q — > P", we must choose the second variant of arrows arrangement. 



6 This differs from the commonly accepted viewpoint that one should instead take the primitive of y(£). But the 
1-differential ydx, being the Whitham differential dS in other terminology [7], [8], plays the fundamental role, not its 
primitive. 
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Figure 1: The arrangement of integration contours on the Riemann surface. 

4 Inverting the loop insertion operator. Free energy 
4.1 The iJ-operator 

We now introduce the operator that is in a sense inverse to loop insertion operator (7). Let 7 

^ J f°°+ n—l - 

H- = -re Soo+ V(x) ■ --re Soo _V(x) ■ -to ■ - ]T S t f -. (47) 

1 1 Joo- ~[ J Bi 

The arrangement of the integration contours see in Fig. 1. We now calculate the action of H on 
the Bergmann bidifferential B(x,q) using again the Riemann bilinear identities. We first note that 
B(x,q) = d x dE X) Q Q {q) and we can evaluate residues at infinities by parts. Then, since dE X) Q {q) is 
regular at infinities, we substitute for V'(x) 2y(x) + 2to/x asu oo + and —2y(x) + 2to/x as x — > oo_ 
thus obtaining 

-resoo + (y(x) + ^ dE XtQo (q)dx + reSoo_ ^-y(x) + ^ dE XtQo (q)dx 

n—l „ 

-t dE X:Qo (q) -E^f B (Q> X )> ( 48 ) 

whence the cancelation of terms containing to is obvious, and it remains to take the combination 
of residues at infinities involving y(x). For this, we cut the surface along A- and £>-cycles taking 
into account the residue at x = q. The boundary integrals on two sides of the cut at Bi then 
differ by dE X: Q (q) — dE x+ r q q {q) = 0, while the integrals on two sides of the cut at Ai differ by 

dE Xt Q (q) — dE x+ r q {q) = §b z B{q, x), and we obtain for the boundary term the expression 

V i y{x)dx i B(q,g), 
S J Ai J Bl 

which exactly cancel the last term in (48). It remains only the contribution from the pole at x = q, 
which is just —y{q)- We have therefore proved that 

H-B(;q) = -y(q)dq. (49) 

Let us now consider the action of H on Wk(-) subsequently evaluating the action of loop insertion 
operator (7) on the result. Note first that the only result of action of d/dV(p) on the operator H 
itself are derivatives dV(x)/dV(p) = — l/(p — x) (and recall that by definition \p\ > |x|, i.e., instead 
of evaluating residues at infinities one should take residues at x = p, and we obtain 

d 



dV(p) 



(H-W k (-)) = W k (p)+H-W k (-,p). (50) 



7 This definition works well when acting on 1-forms regular at infinities. Otherwise (say, in the case of Wo(p)), the 
integral in the third term must be regularized, e.g., by replacing it by the contour integral around the logarithmic cut 
stretched between two infinities. 
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For the second term, due to the symmetry of W k (p, q), we may take the point p as the root of the tree 
subgraph. Then, the operator H always acts on £?(•,£) where £ are integration variables of internal 
vertices. However, if this vertex is an innermost (i.e., there is no arrowed edges coming out of it), then 
the 1-form y(£)d£ arising under the action of H (49) cancels the corresponding form in the integration 
expression, and the residue vanishes being nonsingular at the branching point. If there is an outgoing 
arrowed edge, say d-E Pi Q (£) (can be just one as one line is external), then, again, we can push the 
integration contour for £ through the one for p; the only contribution comes again only from the pole 
at £ = p. The value of the residue is doubled, and we come to the following graphical representation 
for the action of the operator H: 



H- 



H- 



Q 



- Q 



p ; 



Q 



(51) 



For H q -W k {q,p) 



H, 



'i ' dV(q) ^fc(P)> we obtain that for each arrowed edge on which the action of (7) 
he 



produces the new vertex, the inverse action of H q - just give the factor —1 and on each nonarrowed 
edge on which the action of (7) produces the new vertex, the inverse action of H q - just gives zero. As 
the total number of arrowed edges is 2k — 1 for any graph contributing to the sum of diagrams, we 
obtain that 

H q -W k (q,p) = -(2k-l)W k (p) 
and, combining with (50), we just obtain 



° (H q ■ W k {q)) = -(2 - 2k)^J^F k , 



dV(p) 



dV(p) 



(52) 



and, since all the dependence on filling fractions and to is fixed by the claim that the answer depends 
only on [i a and derivatives of M(p) of finite orders at the branching points, we conclude that 



F k = 



2k -2 



HW k . 



(53) 



This is our final answer for the free energy. It permits us to calculate all F k except the contribution 
at k = 1 (torus approximation). It was however calculated by a direct integration in [21]. All other 
orders can be consistently calculated. For this, we only introduce one new vertex o in which we have 

J y(s)ds t 

Although the integral term Jq q y(s)ds is nonlocal, its constant part Jq° drops out 



WO 2ti 



of the residue in the 1MM case, and we can integrate it in the neighborhoods of every branching point 
p a separately. That is, we need only the expansion of this integral as £ — > fi a . Then, say, the genus 
two contribution is provided by the sum of three diagrams 



2F 2 = 2 <£HO + 2 




(54) 



For completeness, we also present 
1 



which can also be written like in (53) 



y(p)y(q) \og\p- q\dpdq 



Fq — —^Hrcg ■ Wq. 



(55) 
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where H reg is explained in many articles (see [18]), and 

Fi = —log ^ M(^)A 4 (deW) 12 ^ , (56) 
where A is the Vandermonde determinant for fi a and a is (n — 1) x (n — l)-matrix 

'« = A, ikt ; 

note that deter = det _1 i?jj, where iJjj are coefficients of the polynomials expressing canonically 
normalized holomorphic 1-differentials: 

A 

dwi = — d£, 

y(0 



and, for instance, 
Meanwhile, it is easy to see that H • W\ is constant. Indeed, 

the first term has simple zero with residue 1/2 to be compensated by the double pole of -£?(£,£) ~ 
4(g_ 1 M p i an< i the total answer is then just the constant n/4. 



4.2 Scaling relation 

We now a little demystify relation (52). Indeed, recall that for any functional F of a finite number of 
"local" variables, which are in our case the branching points fi a and the moments 



(k-l)\ dp 



P=Mq 



we have the relations 



and 



dF 



dF 
dF 



d£ 



dF _ 

dSl ~ f Bi dV(0 
Relation (52) is then equivalent to the relation 



oo . i . n— 1 ., 



an 



F = 



(57) 
(58) 

(59) 



satisfied by the total 1MM free energy. Indeed, coming back to basics of the matrix model, note that 
the operator in (59) just generates the scaling transformations: the part with derivatives w.r.t. t& 
multiplies all vertices by the same scaling factor p simultaneously multiplying propagators by 
The action of derivatives in to and Si results in multiplying all index loops (faces of the fat graph) by 
p. Therefore, for any graph, the total factor is 

# vertices - # edges + # faces 2— 2g 

P Pi 



and it is exactly canceled by the scaling of the formal expansion parameter fi — ► ph. 
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